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Diverse vulnerabilities of Bangladesh’s agricultural sector in 16 sub-regions are assessed using 
experiments designed to investigate climate impact factors in isolation and in combination. Climate 
information from a suite of global climate models (GCMs) is used to drive models assessing the 
agricultural impact of changes in temperature, precipitation, carbon dioxide concentrations, river floods, 
and sea level rise for the 2040-2069 period in comparison to a historical baseline. Using the multi-factor 
impacts analysis framework developed in Yu et al. (2010), this study provides new sub-regional 
vulnerability analyses and quantifies key uncertainties in climate and production. Rice ( aman , boro , and 
aus seasons) and wheat production are simulated in each sub-region using the biophysical Crop 
Environment REsource Synthesis (CERES) models. These simulations are then combined with the MIKE 
BASIN hydrologic model for river floods in the Ganges-Brahmaputra-Meghna (GBM) Basins, and the 
MIKE2 1 Two-Dimensional Estuary Model to determine coastal inundation under conditions of higher mean 
sea level. The impacts of each factor depend on GCM configurations, emissions pathways, sub-regions, and 
particular seasons and crops. Temperature increases generally reduce production across all scenarios. 
Precipitation changes can have either a positive or a negative impact, with a high degree of uncertainty 
across GCMs. Carbon dioxide impacts on crop production are positive and depend on the emissions 
pathway. Increasing river flood areas reduce production in affected sub-regions. Precipitation uncertainties 
from different GCMs and emissions scenarios are reduced when integrated across the large GBM Basins’ 
hydrology. Agriculture in Southern Bangladesh is severely affected by sea level rise even when cyclonic 
surges are not fully considered, with impacts increasing under the higher emissions scenario. 

Published by Elsevier Ltd. 


1. Introduction 

Bangladesh (Fig. 1 ) lies on mostly flat, alluvial land at the mouth 
of the Ganges-Brahmaputra-Meghna (GBM) Basins that drain 
monsoon runoff from a large portion of South Asia, and is widely 
recognized as a country with high sensitivity to climate variability 
and change. Bangladesh uses more than 70% of its land for 
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agricultural purposes (FAOSTAT, 2009), often with multiple 
cropping seasons, and nearly all of the remainder is covered by 
forests, settlements, roads and waterways (MPO, 1986). 
Bangladesh also has high population density, with a current 
population equivalent to half of the population of the United States 
living in an area the size of the state of Iowa. Long-term climate 
threats include a changing distribution of river floods, sea level rise 
in the Bay of Bengal, warming temperatures, and changing rainfall 
patterns, which can exacerbate current vulnerability to climate 
extremes and variations and pose a substantial challenge in 
producing enough food for a growing and developing population 
(Fluq, 2001; Huq et al, 2003; Parry et al., 2004; Asada and 
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Fig. 1. Bangladesh sub-regions as defined for this study, with sentinel weather stations and major rivers. 


Matsumoto, 2009). A beneficial climate impact comes from a key 
driver of global climate change, as elevated carbon dioxide (C0 2 ) 
concentrations enhance photosynthetic production and water 
efficiency in crops to a limited extent (Easterling et al., 2007; 
Fleischer et al., 2010; Kimball, 2010). These processes are among 
the complex challenges facing farmers in Bangladesh, and the 
experiments presented in this study elucidate the agricultural 
responses to projected climate conditions that are unprecedented 
in recent history. 

Bangladesh’s precarious position with respect to climate 
change has elevated awareness of national vulnerabilities both 
domestically and internationally. There have been numerous 
studies of climate change impacts in Bangladesh, but they have 
provided climate projections without quantifying agricultural 
impacts or tended to examine just a subset of the impact factors 
that are inter-compared and combined here. Rahman et al. (2009) 
conducted a comprehensive overview of climate change projec- 
tions for temperature, rainfall, river floods, and sea level rise for 
various Bangladeshi sub-regions and identified adaptation strate- 
gies geared to particular impact factors and vulnerabilities from 
interviews and stakeholder workshops. Similar climate change 
impact reviews are provided by Agrawala et al. (2003), Ahmed 
(2006), and Tanner et al. (2007). Other studies of climate change 
impacts on Bangladeshi rice have focused primarily on coastal 
flooding impacts (Ali, 1999; Sarwar, 2005), temperature and 
carbon dioxide effects (Karim et al., 1994, 1998; Timsina et al., 


1997; Mahmood, 1997, 1998; Timsina and Humphreys, 2006; 
Basak et al., 2009), or river flood impacts (Hassan et al., 2008). 
Much work has also been done to identify and test adaptation 
strategies (Huq et al., 1999; Ministry of Environment and Forests, 
2005; Adger et al., 2007; Thomalla et al., 2005; Ayers and Forsyth, 
2009; Rahman et al., 2009), and in 2009 the Government of 
Bangladesh released its Climate Change Strategy and Adaptation 
Plan (Ministry of Environment and Forests, 2009). 

Government programs (including the construction of embank- 
ments and coastal protection), private enterprise (including the 
expansion of electrified pump-irrigation systems), crop breeders 
(developing drought, flood, or salinity resistant seeds), small farm- 
holders (changing planting dates and/or farm management 
practices) and international partners (e.g. the World Bank, The 
Global Environmental Facility) all have roles to play to build 
efficient climate resilience and adaptive capacity. For a list of past 
and present programs on adaptation in the agricultural sector, see 
Table 7.1 of Yu et al. (2010). Funds are limited, however, so major 
adaptation programs must be carefully prioritized to maximize 
utility. 

This study aims to identify the critical climate factors and 
associated uncertainties in agricultural impacts for each of 16 sub- 
regions in Bangladesh in order to form a basis for prioritizing 
adaptation measures. The impacts of changes in temperature, 
precipitation, C0 2 concentration, river flooding, and sea level rise 
on agricultural production are investigated through simulations 
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where each climate factor is investigated independently. Then, an 
integrated multi-factor analysis allows for sub-regional scale 
examination of the impacts of all climate factors in combination to 
form a framework that may be used to identify and prioritize 
adaptation strategies. This work is an extension of the crop 
modeling portion of a World Bank report on climate change risks in 
the agricultural sector of Bangladesh (Yu et al., 2010), adding 
increased analysis of sub-regional vulnerability and adaptation 
priorities, distinguishing sub-regionally between temperature and 
rainfall impacts, as well as tracking and quantifying sources of 
uncertainty in each sub-region. 

2. Agriculture in Bangladesh 

Agriculture is a pillar of Bangladesh’s economy, using more than 
70% of land area (FAOSTAT, 2009) and accounting for nearly 20% of 
gross domestic product and 65% of the labor force, employed 
primarily on small-holder farms (Yu et al., 2010). Rice production 
occurs on more than 80% of agricultural lands and is grown in three 
growing seasons that span the entire year and are synonymous with 
their rice crop: aman , boro , and aus (BBS, 2005, 2008). Aman grows 
during the monsoon season when rainfall is plentiful. Boro is grown 
during the dry season, after floods have receded, and is restricted to 
irrigated areas. Aus production takes advantage of rainfall during the 
spring transition toward the monsoon that enables a short growing 
season, although the sufficiency of rains varies from year to year. 
Rice cultivation makes up nearly 95% of cereal production in 
Bangladesh, with wheat a prominent dry-season crop grown with 
irrigation (Bangladesh Bureau of Statistics; BBS, 2008). Rice is 
threatened by floods during the pre-monsoon (aus) and monsoon 
(aman) seasons and by heat waves and water scarcity during the dry 
monsoon (boro) season. 


Boro is currently the most productive season, followed by aman, 
with aus considerably lower due mostly to smaller planted area 
(BBS, 2008). Historical cereal production has been increasing for 
decades, rising from about 10 million metric tons (MT) in the early 
1970s to almost 30 MT in 2001 (USDA, 2008). Annual variations 
and long-term trends in agricultural production come from many 
factors, including climate events (particularly floods), resource 
limitations, sociopolitical events, and the implementation of 
modern agricultural practices, so isolating the sensitivity of 
historical agricultural production to climate factors is challenging. 

3. Biophysical approach and model framework 

3.1. Models and methods 

Biophysical process models in this study occupy a range of 
scales, with crop models focusing on representative farms 
throughout Bangladesh, river floods simulated over the entire 
GBM Basins, and coastal inundation determined along the 
Bangladesh coastline. Differences in the computational resources 
available for each model limited the ability to set up identical 
ensembles of climate drivers for each set of model simulations. 
However, the full range of crop model experiments provided a 
framework for integrating the highest possible number of river 
flood and sea level rise simulations and their impact on agricultural 
lands. This modeling framework was described in Yu et al. (2010), 
and is therefore only summarized below and in the Supplementary 
Online Material Appendices. 

Fig. 2 describes temperature and precipitation characteristics of 
the climate scenarios examined in this study, along with baseline 
observations. Climate scenarios are generated according to local 
historic conditions and simulated climate changes from global 
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Fig. 2. Mean monthly, annual, and seasonal (a) temperature (°C), and (b) precipitation (mm/day) for current climate and projected scenarios in the 2050s averaged across all 
16 sub-regions. The starred black line represents the mean of the baseline conditions, while the box and whiskers diagrams represent the 32 GCM/emissions scenario 
combinations for each 30-year period. Note: ‘wet’ season = May-September; ‘dry’ = October- April; ‘aus’ = April-July; ‘aman’ = July-October; ‘boro’ = December-March (this 
is also the wheat season). 
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climate model (GCM) output contributed to the World Climate 
Research Programme’s (WCRP’s) Coupled Model Intercomparison 
Project phase 3 (CMIP3) multi-model dataset (Meehl et al M 2007) at 
the Program for Climate Model Diagnosis and Intercomparison 
(PCMDI; http://www-pcmdi.llnl.gov). These climate simulations 
were analyzed by the studies reviewed in the Intergovernmental 
Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). 
Climate change scenarios are generated by comparing a given GCM’s 
2040-2069 period (referred to as the “2050s”) with control 
simulations of that GCM over a 1970-1999 baseline, and then 
imposing these changes on historical observations (daily rainfall and 
maximum and minimum temperatures from BMD observations; 
solar radiation data and short gaps were filled using the Weather- 
man Utility; Hoogenboom et al., 2003 ). This delta method (described 
more completely in Wilby et al., 2004, and Yu et al., 2010) offsets 
many of the common GCM biases but assumes no change in high- 
frequency variability or the frequency of rain events. Both a 
relatively high (A2) and low (Bl) future emissions pathway were 
analyzed and compared to the baseline period, with crop models 
utilizing representative C0 2 concentrations of 556 ppm, 498 ppm, 
and 345 ppm, respectively, set according to projected or observed 
concentrations in the central year of 30-year the period being 
studied (SRES, 2000). Following Yu et al. (2010), an ensemble of 
climate scenarios were created for each sub-region from 16 GCMs 
and 2 emissions scenarios, capturing a consistent temperature rise 
and wide uncertainty among projected precipitation changes. 

The MIKE BASIN hydrologic model was employed over the 
entire GBM Basins that drain through Bangladesh to simulate river 
floods for each baseline and future year. Flood protection 
infrastructure in Bangladesh was also taken into account by a 
Bangladeshi flood model (Nishat and Rahman, 2009; Flopson and 
Webster, 2009). Coastal inundation in the 2050s was simulated by 
the UK Department for Environment, Food, and Rural Affairs in 
collaboration with IWM and CEGIS (DEFRA, 2007) using the 
MIKE21 Two-Dimensional Estuary Model and mean sea levels of 
27 cm (~A2) and 8 cm (~B1) above baseline. More details about 
the flood models are provided in Appendix II of the Supplementary 
Online Material. 

Process-based crop model simulations were run with the Crop 
Environment REsource Synthesis (CERES) rice and wheat models 
from the Decision Support System for Agrotechnology Transfer 
(DSSAT V4.5.0.030; Eloogenboom et al., 2003; Jones et al., 2003). 
CERES models require detailed soil profiles, farm-level manage- 
ment practices, genetic coefficients describing the cultivar, and 
daily meteorological conditions (rainfall, solar radiation, and 
maximum and minimum temperature). CERES models have been 
applied extensively in Bangladesh, although less frequently for 


climate change applications (e.g., Karim et al., 1996; Karim et al., 
1994, 1998; Elussain, 2006; Basak et al., 2009). These models 
include the beneficial effects of enhanced carbon dioxide 
concentrations on plant growth and are sensitive to the 
phenological stage of plant development when particular stresses 
occur, but the exact relationships remain uncertain. For this study 
the CERES model was calibrated for Bangladeshi sub-regions using 
soil and cultivar information provided by S.G. Elussain (2008), 
additional soil profiles drawn from Brammer (1996), and farm- 
level management according to practices recommended by the 
Bangladesh Rice Research Institute (BRRI, 2007). 

Crop model simulations for this study utilize modern cultivars 
and relatively high fertilizer application rates and assume pest-, 
salinity-, and disease-free conditions. This management configu- 
ration follows the recommendations of the Bangladesh Rice 
Research Institute (BRRI, 2007; presented in Appendix I in the 
Supplementary Online Material), but is a more intensive manage- 
ment than is typical for many farms in the region. Mahmood et al. 
(2003) noted the large gap between actual and potential yields in 
Bangladesh, with actual yield averages approximately 1 /6th of the 
potential yields produced under unstressed, high-input conditions 
protected from floods. Projections of future yields provide an 
indication of the magnitude and direction of climate impacts 
which will occur on top of the anticipated reduction in this yield 
gap through economic development. 

4. Definition of sub-regions 

Model and observational data for this project’s impacts assess- 
ments employ different spatial scales and coverage densities. To 
overcome the resulting challenge in consistently integrating 
simulations of different climate impact factors, a common spatial 
scale was created using available data and resources that groups 
similar farming areas into 16 sub-regions that cover all of 
Bangladesh. These sub-regions form the common spatial scale on 
which climate change factors are compared. Sub-regions were 
defined according to climatology, agricultural use, terrain, and 
political boundaries, making them very similar to agro-ecological 
zones with the added consideration of introducing political 
boundaries to facilitate clear communication and appropriate 
jurisdiction over adaptation strategies put into action. These sub- 
regions are similar to those developed by Mahmood et al. (2004), but 
were chosen so as to insure that each had an observation station of 
the Bangladesh Meteorological Department (BMD) with observa- 
tions ideally going back to 1970. The number of sub-regions was 
limited to be consistent with simulation resources. Fig. 1 shows the 
16 sub-regions as defined by the Center for Environment and 


Table 1 

Sub-regional climate, soil, planted area, and flood-protected area information, as well as the sub-regions included in the simulations of the Basins (bold) and coastal 
(italicized) flood models. 


Sub-region 

BMD climate station 

Soil profile location 

Aus area (ha) 

Aman area (ha) 

Boro area (ha) 

Wheat area (ha) 

Flood-protected area (%) 

1 

Dinajpur 

Dinajpur 

65401 

613853 

543438 

174969 

45 

2 

Rangpur 

Rangpur 

42491 

815334 

976956 

149176 

48 

3 

Ishwardi 

Jessore 

21165 

116985 

86369 

53228 

76 

4 

Tangail 

Karatia 

46826 

239834 

261113 

48131 

15 

5 

Dhaka 

Ghatail 

13458 

51661 

32107 

8225 

9 

6 

Mymensingh 

Phulpur 

100018 

610184 

561919 

39339 

17 

7 

Sylhet 

Biani Bazar 

124288 

350173 

589900 

10034 

27 

8 

Srimangal 

Srimangal 

24304 

64558 

19294 

541 

22 

9 

Comilla 

Shalpur 

129498 

330897 

426250 

33260 

45 

10 

Chittagong 

Chittagong 

38209 

23194 

74850 

39 

31 

11 

Rangamati 

Srimangal 

30320 

14628 

71881 

17 

2 

12 

Maijdee Court 

Hatiya 

45945 

121579 

26894 

103 

58 

13 

Jessore 

Jessore 

98834 

359053 

333375 

67317 

41 

14 

Faridpur 

Jessore 

167592 

351018 

286413 

49575 

38 

15 

Patuakhali 

Satkhira 

240201 

198407 

162463 

6399 

48 

16 

Khulna 

Satkhira 

14809 

215747 

78225 

1794 

35 
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Table 2 

Climate information for each sub-region. The representative BMD station, its code, and annual mean climate statistics during the 1970-1999 baseline period. 


Sub-region 

BMD station location 

BMD station Code 

Mean T max ( C) 

Mean T min (°C) 

Mean rainfall (mm) 

Mean sunshine (MJ/m 2 /day) 

1 

Dinajpur 

10120 

30.1 

19.7 

2003 

16.9 

2 

Rangpur 

10208 

29.7 

19.9 

2239 

17.2 

3 

Ishwardi 

10910 

31.0 

20.3 

1652 

17.3 

4 

Tangail 

41909 

30.3 

20.8 

1902 

16.7 

5 

Dhaka 

11111 

30.6 

21.6 

2148 

17.6 

6 

Mymensingh 

10609 

30.0 

20.5 

2255 

16.4 

7 

Sylhet 

10705 

29.6 

20.2 

4150 

16.7 

8 

Srimangal 

10724 

30.4 

19.4 

2421 

17.0 

9 

Comilla 

11313 

30.1 

20.9 

2054 

17.2 

10 

Chittagong 

11921 

30.2 

21.6 

2931 

18.0 

11 

Rangamati 

12007 

30.2 

21.4 

2532 

17.0 

12 

Maijdee Court 

11809 

29.8 

21.6 

3103 

16.9 

13 

Jessore 

11407 

31.4 

20.9 

1600 

17.2 

14 

Faridpur 

11505 

30.4 

21.1 

1967 

17.2 

15 

Patuakhali 

12103 

30.3 

21.9 

2704 

15.6 

16 

Khulna 

11604 

31.1 

21.6 

1812 

17.4 


Geographic Information Services (CEGIS) in conjunction with the 
rest of the study team. 

Agricultural information about each sub-region, including the 
area growing aus, aman, boro, and wheat, is presented in Table 1. 
Simulated production in each sub-region is calculated by 
multiplying the simulated per hectare yield by the agricultural 
area in each sub-region producing the crop of interest. Agricultural 
areas were assumed to be constant over the coming century as very 
little arable land is not currently farmed, and areas are assumed to 
be evenly distributed across any given sub-region because of the 
lack of finer-scale information. Mean baseline period climate 
variables for each sub-region are presented in Table 2. 


5. Climate factors 

5.1. Direct climate factors 

Direct climate factors are those experienced directly at the farm 
level from the atmosphere, including changes of temperature, 
precipitation, C0 2 concentrations, and solar radiation that affect 
crop development. 

5.1.1. Temperature 

Warmer conditions associated with climate change lead to a 
higher rate of respiration in growing crops, requiring greater 
amounts of carbon for consumption, speeding up phenological 
development, and increasing the rate of tissue aging. Temperature 
changes also affect the timing of seasons, increase potential 
evapotranspiration, and may push plants closer to damaging 
thresholds (Solomon et al., 2007; Easterling et al., 2007; Hatfield 
et al., 2008). The ability of a crop to withstand heat stress depends 
upon the developmental stage when a heat wave occurs, which is a 
focus of ongoing research and a sensitivity not completely 
captured by the CERES models. 

To investigate the effects of isolated temperature changes, AT 
climate scenarios were generated for the CERES crop models to be 
equal (in radiation, precipitation, and C0 2 ) to the 1970-1999 
baseline period observations, but with daily maximum and 
minimum temperatures modified according to mean local monthly 
changes in each GCM’s simulated temperature (i.e., the delta 
method applied for temperature only; Wilby et al., 2004). For these 
scenarios the baseline temperature was modified according to the 
difference between model means for the 2040-2069 period and 
the means for the 1970-1999 baseline period for each GCM. AT 
climate scenarios were created for all 16 GCMs for both the A2 and 
B1 emissions scenarios. At this stage of the analysis, no flood 
damages were applied. 


5. 1 .2. Precipitation 

Aus and aman are rainfed crops that periodically experience 
drought stress (particularly in the northwest). Precipitation 
changes in future climate simulations vary widely across models 
and regions, making them among the most uncertain climate 
change projections. AP climate scenarios isolating the effect of 
precipitation changes were created for the CERES crop models to be 
equal (in radiation, minimum temperature, maximum tempera- 
ture, and C0 2 ) to the 1970-1999 baseline period observations, but 
with daily precipitation totals modified according to mean local 
monthly percentage changes in each GCM’s simulated precipita- 
tion (i.e., the delta method for rainfall only; Wilby et al., 2004). 
These percentage changes are the difference between the GCM’s 
2040-2069 period values and those for the 1970-1999 baseline 
period for that GCM. AP climate scenarios were created for all 16 
GCMs for both the A2 and B1 emissions scenarios, and no flood 
damages were applied. Because the boro and wheat cropping areas 
are assumed to have unlimited irrigation, the impacts of 
precipitation changes are negligible. 

5.1.3. C0 2 

Biophysical crop models, controlled chamber experiments, and 
Free Air C0 2 Enrichment experiments (FACE; Hendry and Kimball, 
1 994) have demonstrated the relevance of increased C0 2 concen- 
trations on agricultural production, but the extent to which large- 
scale field crops will respond to C0 2 is the subject of considerable 
debate (Easterling et al., 2007; Long et al., 2006; Tubiello et al., 
2007a, b; Ainsworth etal.,2008; Hatfield etal.,2008; Kimball, 2010; 
Fleischer et al., 201 0) and ongoing research (Boote et al., 201 0). C0 2 is 
a primary element of photosynthesis, and plants respond to elevated 
levels by increasing the rate of primary production. High C0 2 
concentrations also increase root densities, allow a plant to make 
more efficient gaseous transfers with its environment, and increase 
water efficiency (Hatfield et al., 2008). Additional research is needed 
to determine the agricultural effects of C0 2 concentrations for more 
diverse locations, a wider variety of cultivars, and across fields with a 
range of nitrogen applications. 

To isolate the impact of enhanced C0 2 concentrations, AC0 2 
climate scenarios were generated for the CERES crop models 
that were identical to the 1970-1999 baseline observations 
(radiation, minimum temperature, maximum temperature, and 
precipitation), but where 1985 C0 2 concentrations (345 ppm; 
Keeling and Whorf, 2005) were replaced by the center (2055) 
value of the 2040-2069 period. Because all of the GCMs 
simulating a given emissions scenario use the same C0 2 
concentrations, AC0 2 climate scenarios were created according 
to emissions scenario (A2 = 556 ppm and B1 = 498 ppm) rather 
than for each GCM. 
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Fig. 3. Maximum flood extent (as percentage of total land area) in Bangladesh from 1954-2007 (Rahman et al., 2009). Note that missing years may not have been flood-free, 
and that Bangladesh’s total area is ~ 144,000 km 2 . 


5.2. Indirect climate factors 

Crops may also be affected by climate events that originate 
away from the farm but result in local flooding. 

5.2.1. GBM river floods 

The monsoonal discharge of water from the GBM Basins results 
in regular seasonal floods affecting Bangladesh with extents 
varying by year and in response to particular storm events. 
Inundation of rice crops suffocates gas exchanges necessary for 
photosynthesis and respiration and can also uproot plants or wash 
away top soil, reducing overall yields. Fig. 3 shows the areal extent 
of national floods over recent decades, showing large variations in 
total extent, and also events that covered nearly 70% of the total 
land area (in 1998; Rahman et al., 2009). 

Flood damages were calculated using a method adapted from 
S.G. Hussain, with damages depending on the duration of the flood, 
the fraction of the plant inundated, and the developmental stage of 
the crop when the flood strikes (Hussain, 1995; described in 
Appendix II of the Supplementary Online Material and hereinafter 
referred to as the Hussain method). Mean plant heights and 
phenological stages for each calendar day were taken from CERES- 
Rice simulations calibrated for the 1970-1999 baseline period 
under observed meteorological conditions. All land was classified 
for flood depth, taking areas protected by polders and other flood 
embankments into account. Simulations were run for the 1 1 sub- 
regions entirely in the GBM Basins (bold names in Table 1) using 
the MIKE BASIN hydrologic model at the Institute for Water 
Modeling (Nishat and Rahman, 2009). To isolate the effects of 
changing GBM river flood extents, AFlood climate scenarios were 
generated that used baseline values for the direct variables 
(temperature, precipitation, radiation, and C0 2 ), but used modified 
climate in the GBM Basins for input to the MIKE BASIN hydrologic 
model. The catchment climates driving river flood simulations 
were based on 1978-2007, which was the closest 30-year period 
where GBM conditions were available. These were modified 
according to each GCM’s mean monthly change in simulated 
temperature and percentage change in simulated precipitation 
(2040-2069 period in comparison to the 1970-1999 baseline 
period from that GCM). Differences between the 1970-1999 
baseline period used for the crop model scenarios and the 1978- 
2007 baseline period used for the river flood models were small in 
comparison to the differences between 2040-2069 and either of 
the baseline periods. Due to resource constraints, scenarios from 
only 5 GCMs (GFDL CM2.1, Delworth et al., 2006; MIROC3.2, 
medium resolution, Hasumi and Emori, 2004; MPI Echam5, 
Jungclaus et al., 2006; NCAR CCSM3.0, Collins et al., 2006; and 
HadCM3, Johns et al., 2006) were developed for the A2, Bl, and 
20th century (for baseline conditions) scenarios during the 
monsoon seasons; to our knowledge this is a larger ensemble 
than has been previously produced for Bangladesh. No sea level 
rise or direct climate impacts were simulated for crop production 


at this stage of the analysis. Additional details on the river flood 
models, the flood damage projections, and the adaptation of the 
Hussain method for multi-factor model projections are provided in 
Appendix II of this article’s Supplementary Online Material. 

5.2.2. Sea level rise 

Sea levels rise in response to thermal expansion of the oceans, 
the addition of meltwater from polar caps and alpine glaciers, and 
local effects including circulation changes and subsidence, erosion, 
and accretion of land. The resulting net encroachment onto 
agricultural land is exacerbated by tides and storm surges, which 
push damaging waves further inland and increase soil salinity. 
Thus, land areas that were inundated by at least 30 cm of seawater 
according to the IWM/CEGIS Bay of Bengal Model (DEFRA, 2007) 
were assumed not to be available for cereal production. Scenarios 
were created for the Bay of Bengal model simulations with mean 
sea level increased to mean 2050 estimates from the IPCC Third 
Assessment Report (Houghton et al., 2001 ; McCarthy et al., 2001), 
and inundation was modeled in 8 coastal sub-regions (italicized 
names in Table 1 ). Approximating across coupled GCM projections, 
DEFRA (2007) created two scenarios for sea level rise above 2000 
levels (A2 = +27 cm; Bl = +8 cm). Global mean sea level in 2000 
was ~3.5 cm higher than during the 1 970-1999 baseline period, so 
these scenarios slightly underestimate the changes from the 
baseline period used by the crop model experiments (Church and 
White, 2006). To isolate the effects of these sea level rises, ASLR 
scenarios were applied to crop model simulations using the 
baseline record of direct climate variables (temperature, precipi- 
tation, radiation, and C0 2 ) without GBM river flooding. Additional 
details about sea level rise damage calculations are described in 
Appendix II in the Supplementary Online Material, including a 
discussion of uncertainty and the likelihood that these projections 
are optimistic due to an omission of cyclone storm surge (Karim 
and Mimura, 2008). 

6. Individual climate change factor impacts 

6.1. Direct climate factor scenarios 

Fig. 4 presents the results of the AT climate scenario 
experiments, showing reductions in yield from warmer tempera- 
tures in all seasons and nearly all sub-regions. The largest declines 
are projected for the irrigated boro and wheat crops. Wheat 
experiences the largest median declines from warmer tempera- 
tures (-12.4% nationally), with declines of more than 15% in sub- 
regions SR-1 (the largest at -19%), SR-6, SR-11, and SR-15. SR-1 is 
the most productive sub-region for wheat, while the latter two 
produce very little wheat at present. Among the ensemble of crop 
model simulations driven by scenarios from 1 6 GCMs, median boro 
production decreases by 12.1% nationally, with losses of greater 
than 1 0% in 1 1 sub-regions, including 6 of the top 7 boro-producing 
sub-regions. Aman is projected to be the crop least affected by 
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b) Aman 


a) Aus 


d) Wheat 





Fig. 4. Median sub-regional production impacts for the AT experiment for (a) aus, Fig. 6. Median sub-regional production impacts for the AC0 2 experiment for (a) aus, 
(b) aman, (c) boro, and (d) wheat. (b) aman, (c) boro, and (d) wheat. 


temperature (median -1.5% national production change), particu- 
larly in the highly productive northwest (SR-1, SR-2) where 
temperatures are slightly cooler than in most of Bangladesh. 
Median projections show a -7.3% median change in national aus 
production, distributed fairly evenly across nearly all sub-regions, 
although temperature changes in SR-7 cause more than a 12% loss 
in production. 

The results of the A P climate scenario experiments for rainfed 
aus and aman are presented in Fig. 5. Projected rainfall changes 
reduce median national aus production by 1.4%, with losses in all 
sub-regions except SR-7 (which shows slight gains of 1.5%). The 
largest declines occur in SR-4 (nearly 5%), although the most 
productive aus region (SR-1 5) declines by only 1.2%. Aman 
simulations were barely sensitive to the projected changes in 



Change in production (%) 

Fig. 5. Median sub-regional production impacts for the A P experiment for (a) aus, 
(b) aman. Note: irrigated seasons are not shown. 


rainfall, with only SR-3 and SR-4 near a 1% median change in 
production. In general, rainfed aus and aman are not strongly 
limited by the mean seasonal rainfall, but are likely to be more 
sensitive to shifts in the variability of damaging dry spells and 
delayed monsoon onset. These changes, as well as changes in the 
demand and availability of groundwater (from pumping and 
natural fluxes) for boro and wheat irrigation, merit further study 
(e.g., Shahid, 2011; Arnell et al., 2011). 

AC0 2 scenario results are presented in Fig. 6. As expected, all 
seasons and sub-regions show gains in production under higher 
C0 2 concentrations with median national production increases of 
11.7% (aus), 4.9% (aman), 8.3% (boro), and 10.1% (wheat). Sub- 
regional responses vary, as C0 2 concentration changes interact 
with the soil and weather patterns for sub-regional heterogeneity. 
Probably reflecting an added benefit from increased water 
efficiency in plants in higher C0 2 concentrations, rainfed aus 
production gains are the largest among rice seasons. Every sub- 
region sees at least an 8% increase in aus production; SR-1 5, which 
has the largest aus area, leads with more than 15% gains. In 
contrast, the wet aman season does not often suffer from water 
stress and loses more nitrogen to leaching which makes it less able 
to take full advantage of C0 2 enhancement, with gains between 
2.5-7% in all sub-regions. The irrigated boro and wheat seasons 
have less leaching and show large gains from enhanced C0 2 , 
indicating that primary production effects are large. SR-6, SR-1 5, 
and SR-1 6, which were the most sensitive to C0 2 concentration 
changes, had the highest levels of NH 4 in their soil profiles, 
indicating additional nutrient limitations to C0 2 effects. 

6.2. Flood hazards 

The effect of mean climate changes on river flood damages are 
presented in Fig. 7a and b for aus and aman (the dry season crops 
were assumed to be flood-free). Simulated river flood damages 
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Fig. 7. Median sub-regional production impacts for the AFlood experiment for (a) 
aus, and (b) aman. Note: Boro and wheat were assumed to be flood-free, and sub- 
regions in white were not simulated. 


only reduce national aus production by a median of 1.9%, although 
sub-regions in northern Bangladesh bear the brunt of the early 
season floods, with the worst damage to the east (Fig. 7a). The 
highest median increase in flood damages occurs in SR-7 (8.9% loss 
in production), where low-lying areas and flash floods are 
common. 

During the later aman season (Fig. 7b), substantial climate 
change increases in flood damage occur in nearly all sub-regions 
simulated by the hydrologic model, lowering national production 
by a median of 3.2%. Increases in damage are generally greater 
towards the east, with SR-7 (9.2% loss) and SR-4 (8.2% loss) 
showing the largest median damages; however, the eastern hills in 
SR-8 are least affected. Although SR-9 was not simulated by the 
MIKE BASIN model, its banks of the Meghna are also likely to 
experience substantial flood damage, similar to the affected sub- 
regions that surround it. The mean impact of changing river floods 
from climate change is much smaller than typical year-to-year 
flood variations (Fig. 3), but climate change is projected to raise the 
level of the average flood and make it more likely that extreme 
events will be damaging. Climate change-induced effects on the 
frequency and severity of extreme floods merit further study, as 
more extreme interannual variability would dwarf the river flood 
impacts of mean climate shifts alone (Mirza, 2002; Yu et al., 2010; 
Mirza, 2011). 

Fig. 8 shows projections of production lost to coastal 
inundation associated with mean sea level rise (average of A2 
and B1 scenarios). The same percentage of land is lost from the 
coastal sub-regions in each season, but the dry season crops are 
disproportionately grown away from the coast, so they are 
affected less. National losses of 2.3% (aus), 1.3% (aman), 0.7% 
(boro), and 0.2% (wheat) are projected when the A2 and B1 
scenarios are averaged. Production in SR-16 (-13.3%) and SR-1 5 (- 
8.7%) is most affected. As the current world emissions trajectory is 
higher even than the A2 scenario and sea levels are rising more 
rapidly than previously thought due to complex ice sheet 
dynamics (Alley et al., 2005), it is useful to note that the A2 
scenario simulations project much higher coastal damages (SR- 
16 =-20%; SR-15 = —13.1%; SR-10 = —10%). Simulations of a 
62 cm rise in mean sea level (estimated for the A2 2080s; DEFRA, 
2007) project damages to production because of area loss in 
excess of 31% in SR-15 and nearly 40% in SR-16. With additional 
vulnerability to severe coastal storms, land subsidence from 
groundwater extraction and salinity intrusion, these results 
underscore the long-term threat of sea level rise in Bangladesh 
(Flassan et al., 2006). 
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Fig. 8. Mean sub-regional production impacts for the ASLR experiment (average of 
A2 and B1 scenarios). Note: Sea level rise yield impacts were assumed to affect all 
seasons equally, and white sub-regions were not simulated. 

7. Integrated sub-regional impacts 

7.1. Median impacts 

Multi-factor analyses of production changes were conducted 
for each crop and sub-region, with new scenarios combining 
changes in the direct climate factors (future temperature, rainfall, 
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Fig. 9. Median sub-regional production impacts considering all climate factors for 
(a) aus, (b) aman, (c) boro, and (d) wheat. 
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Table 3 

Uncertainty analyses for multi-factor production impacts. For each crop and sub-region, the median change across all 32 GCM/emissions scenario simulations, the absolute 
difference in the median between the A2 and B1 emissions scenarios (ASRES), and the average inter-quartile range among GCMs in each emissions scenario (AGCM) are 
presented. Shaded boxes indicate that this uncertainty metric makes the sign of the median change difficult to distinguish from zero. 



Aus Rice 

Aman Rice 

Boro Rice 

Wheat 

Sub- 

region 

Median 

Change 

ASRES 

AGCM 

Median 

Change 

ASRES 

AGCM 

Median 

Change 

ASRES 

AGCM 

Median 

Change 

ASRES 

AGCM 

1 

1.5 

0.5 

5.3 

3.5 

0.8 

1.2 

-5.4 

0.8 

8 

4 

1.7 

7.8 

2 

-0.7 

2.3 

7 

2.7 

0.5 

2.1 

-6.8 

0.7 

6.2 

3.9 

1.3 

4.1 

3 

2.9 

0.8 

4.3 

-0.8 

1 

2.5 

-2.8 

1.4 

8.3 

4.8 

1.7 

2.5 

4 

-0.6 

3.4 

6.8 

-4.6 

1.8 

4.8 

-1.5 

1.1 

7.2 

1.9 

0.4 

2.7 

5 

-0.1 

1.5 

9.2 

0.5 

0.3 

0.9 

-4.8 

0.3 

3.4 

5.2 

1.8 

6 

6 

0.7 

0.9 

5.7 

-2.4 

1.3 

2.2 

-4.3 

1.5 

6.2 

7.6 

3.2 

8.1 

7 

-11 

7 

5.4 

-8.3 

6.4 

1.1 

4.8 

0.7 

2.9 

4.1 

2.1 

3.1 

8 

-1.1 

1.1 

3.9 

2.7 

0.5 

1.6 

-9.2 

1.1 

10.3 

4.7 

2.3 

3.4 

9 

-0.3 

0.7 

3.8 

3.4 

0.8 

2.8 

-7.1 

0.9 

6.1 

-0.9 

1.6 

3.1 

10 

-1.2 

1.7 

4.7 

1.7 

2.3 

1.1 

-7.4 

3.4 

2.3 

0.6 

2.1 

3.3 

11 

1.7 

1 

3.3 

3.5 

0.6 

1.4 

-5.4 

0.3 

3.4 

5.1 

4.1 

6.1 

12 

-5.1 

5.7 

3.7 

-0.7 

4.4 

1.2 

-10.5 

5.4 

3.8 

-6.7 

6.6 

2.8 

13 

3.4 

1.5 

3.6 

2.8 

3.8 

5.9 

-5.8 

1 

5.4 

2.2 

0.5 

2.7 

14 

2.4 

2 

4 

-1.8 

2.9 

5.9 

-6 

0.3 

5.2 

3.6 

0.9 

2.8 

15 

-3 

7.5 

5.8 

-9.6 

7.6 

0.5 

-12.1 

9.8 

5.9 

1 

3.6 

9.1 

16 

-10.5 

13.9 

4.1 

-10.8 

12.8 

1.3 

-20.8 

12.9 

4.5 

-12.6 

12.6 

6.3 


and C0 2 concentration changes) with river flood changes and 
coastal inundation from sea level rise applied as multiplicative 
factors (Fig. 9 and Table 3). As expected, the effects of simultaneous 
changes in rainfall, temperature, and C0 2 concentrations are not 
the sum or product of each individual factor’s effect (Figs. 4-8), as 
biophysical interactions between direct climate variables are often 
non-linear. For example, elevated C0 2 concentrations can partially 
mitigate the impact of increased drought, while a combined 
drought and heat wave can cause more damage than either factor 
on its own. Sea level rise has a substantial, year-round impact on 
coastal regions, with damages in SR-16 strongest for each crop. 
Fig. 8 shows acute sub-regional contrasts that may be overlooked if 
only the national production is considered. 

Aus, in particular, shows a wide contrast between sub-regions 
with substantial gains (SR-3, SR-13, SR-14) and those with severe 
losses (SR-7, SR-16; Fig. 9a). The gaining regions are not strongly 
affected by increases in early-season river floods, and have a 
favorable balance of potent C0 2 enhancement and moderate stresses 
from increased temperatures. Aus declines in SR-1 6 are due primarily 
to sea level rise; SR-7 is projected to suffer from the most severe 
temperature and flood impacts. The historically most productive aus 
region (SR-1 5; Table 1) only experiences a median loss of 3% due to 
the strongest C0 2 aus enhancement among the sub-regions. 

The simulated sub-regional pattern of aman changes reflects 
the dominance of the indirect climate change factors; river 
flooding and coastal inundation (Fig. 9b). The effects of the direct 
climate change factors, on the other hand, are relatively less 
important and are consistent across all sub-regions. The largest 
aman damages occur in the coastal region (SR-1 5, SR-16), 
followed by often-flooded SR-7 and SR-4 where waters from 
the Ganges and Brahmaputra meet. Positive changes in SR-9 are 
probably optimistic because flood damages were not considered 
there. Aman production increases in the historically most 
productive sub-region (SR-1), but projections of a gain in SR-2 
are somewhat offset by a loss in SR-3 (the next two most 
productive sub-regions). 


Boro simulations show substantial median declines in all sub- 
regions except SR-7, which experiences the weakest temperature 
impact and a strong C0 2 enhancement in the part of the country 
with the most rain and the lowest average maximum temperatures 
(Fig. 9c). Production losses are the largest on the coast; however, 
losses are nearly 7% in SR-2 in the northwest, which is by far the 
most productive sub-region for boro. 

Wheat production is projected to increase in nearly all sub- 
regions with substantial historical production (Fig. 9d), although 
sea level rise affects coastal areas. ~4% gains are projected for the 
two most productive wheat sub-regions, SR-1, SR-2 in the 
northwest; SR-1 experiences strong C0 2 and temperature effects 
while both factors are slightly moderated in SR-2. SR-9 is the only 
adversely affected sub-region that currently grows large amounts 
of wheat, owing in part to coastal inundation from sea level rise. 
Changes in rainfall and river floods have little effect on boro and 
wheat growth (as they are irrigated in the dry season), but 
hydrologic effects on groundwater recharge during the monsoon 
season could have a lasting impact on irrigation water availability 
that merits further study in a climate change context. 

7.2. Implications of uncertainty 

The integrated multi-factor modeling framework described 
above enables a detailed examination of climate impacts uncertain- 
ty in the impacts assessment process for Bangladesh. It must be 
emphasized, however, that the uncertainties described reflect the 
simulation process rather than the true overall predictive uncer- 
tainty for Bangladesh, which is also dependent on many factors 
(known and unknown) not considered in this study. 

Table 3 shows two primary sources of framework uncertainty in 
the combined impacts for each sub-region. The first is related to 
the differences between societal pathways (via the A2 or B1 
emissions scenarios), and the second is related to the differences 
between the many GCMs simulating future climate. The former 
depends on policy makers deciding and enforcing international 
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emissions targets, while the latter indicates an imperfect model 
representation that may be improved through continuing GCM 
development. For each sub-region and crop the following are 
shown: (1) the median of the 32 futures (16 GCMs x 2 SRES 
emissions scenarios); (2) ASRES, or the difference between the 
median of the 1 6 A2 simulations and the 16 B1 simulations; and (3) 
AGCM, or the average range of the middle 50% of the 16 GCM 
simulations for each emissions scenario. Thus, this table allows 
projected changes to be put in the context of uncertainty, although 
the relatively small sample size restricts the application of more 
robust statistical tests. 

In general, uncertainty among GCM projections is larger than that 
from emissions scenarios, although of course the metrics are 
different and therefore only broad comparisons may be made. 
ASRES is largest in coastal regions across all crops, reflecting the 
importance of societal pathway in the severity of sea level rise 
damages. Although direct climate factors are not strongly sensitive 
to emissions scenarios in this time period, the extent of sea level rise 
depends heavily on changes in the more sensitive Arctic and 
Antarctic (Shepherd and Wingham, 2007). In SR-7 aus and aman are 
also highly sensitive to the emissions scenario, which affects the 
severity of river floods on the Meghna. For all factors (both those 
with positive and negative production implications), A2 changes are 
further from zero than the B1 changes due to larger deviations in 
greenhouse forcing. AGCM is dependent on the cropping season, 
with smaller values (on average) for aman, which is more robust to 
projected changes in direct climate factors, than for other crops. 

Many sub-regions’ median projected production changes are not 
statistically separated from zero, as they are closer to zero than half 
of either ASRES or AGCM. This reflects a substantial number of 
projected outcomes where the change is of the opposite sign (shaded 
in Table 3). This is particularly striking for aus, where the strength 
and timing of pre-monsoon rains and the monsoon onset varies 
considerably among GCMs, affecting temperatures, rainfall, and 
flood damages (Bollasina and Nigam, 2009). Aman impacts are quite 
robust despite relatively small median changes, while boro and 
wheat changes are mostly large enough to be clearly distinguished 
from zero. For all sub-regions and crops, Table 3 also demonstrates 
that the range of simulations includes projections that are 
considerably more severe than the median change. As the 
vulnerability of many communities increases non-linearly with 
the scale of an impact, the full range of projected changes should be 
examined in designing and prioritizing adaptation strategies. 

The individual climate factor impacts in Section 3 may also be 
used to analyze the sensitivities of combined impacts to errors in 
each factor. For example, as research continues to reveal the extent 
to which C0 2 is beneficial to tropical rice and wheat production 
(Hatfield et al., 2008), small changes are likely to have profound 
implications on these projections as C0 2 is the only beneficial 
factor for most sub-regions and crops. Likewise, insights suggest- 
ing increased variability in river floods (e.g. Mirza, 2002), which are 
not modeled here, could elevate the importance of that climate 
change factor. Thus, the present work provides a good basis for 
continuing evaluation of the expected consequences of climate 
change on agriculture in Bangladesh. 


8. Priority targets for adaptation 

Projections of individual and integrated climate change impact 
factors reveal the need for adaptation plans on seasonal and sub- 
regional (rather than annual and national) levels. These projections 
aid in assessing the potential benefit of an adaptation strategy, but 
the costs of implementation and maintenance will also be major 
considerations (see Yu et al., 2010, for an extensive list of potential 
adaptation strategies). 



\ 


Fig. 10. Regions of Bangladesh with common climate change factors of importance. 
Farms in these regions face similar challenges from climate change and may be 
particularly interested in sharing adaptation strategies. Note that the hill tracts sub- 
regions (SR-11 and SR-8) do not fit these larger regions. 

Adaptations increasing climate resilience in the boro season 
should target the highest vulnerabilities as stresses from warmer 
temperatures and limited irrigation are likely to affect this most 
productive crop across the country. Even accounting for projected 
negative stresses from climate change impacts, areas introducing 
irrigation for boro production would increase total production and 
the ability to offset particularly damaging aman flood losses. Aman 
appears to be the crop most resilient to direct climate changes, so 
damages targeting flood protection are projected to be the most 
beneficial for the aman season. Aus projections are the most 
uncertain, but the concentration of aus production in SR-15 is a 
major vulnerability that might suggest some coastal protection. 
Should irrigation water be available, expanded wheat production is 
projected to benefit greatly from elevated C0 2 concentrations in the 
first half of this century. 

In general, three regions of Bangladesh emerge as having 
commonalities in the most impactful climate change factors 
(coastal, river flood, and heat-stressed; Fig. 10). The coastal region 
is most threatened by sea level rise, and a complex and changing 
coastline makes permanent protection infeasible. Salinity intru- 
sion and coastal storms add to the mean sea level rise impacts 
modeled here. The impact of climate change on the river flood 
region is most dependent on the ways in which monsoon rain and 
Himalayan snowmelt flow through Bangladesh, with changes in 
the frequency and intensity of large floods adding to the mean 
flood changes simulated in this study. Polders and other flood 
protections could target this vulnerability. In the heat-stressed 
region, sea level rise and changes in river floods are less impactful, 
so regional efforts to improve crop resilience to warm tempera- 
tures are projected to show the largest local benefit. 

9. Implications and continuing work 

This study has focused on understanding and targeting new 
vulnerabilities introduced by global climate change in Bangladesh. 
Benefits may also be realized through adaptations that increase 
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resilience to current climate variability, which is likely to obscure 
the longer-term climate change signal for the next several decades. 

By the 2050s, climate change is projected to reduce production 
in Bangladesh for all three rice seasons, with only wheat (<5% of 
current national cereal production) showing gains. The boro 
season is expected to be most severely impacted (even when 
irrigation is unconstrained for existing boro areas). C0 2 enhance- 
ment, the only climate change factor with positive effects, will 
likely show diminishing benefits further into the future (Easterling 
et al., 2007), suggesting that detrimental climate change impacts 
will accelerate. Year-to-year climate variability, which exists on 
top of the more gradual changes due to greenhouse gas emissions, 
will continue to produce variations in agricultural production (for 
example between normal years and those with heavy floods or 
strong land-falling cyclones) that are larger than the climate 
change signal. Only by analyzing a longer period of production (as 
in the 30-year series examined here) will the climate change signal 
emerge. 

Over the coming decades many other factors will affect 
Bangladesh’s cereal production and obscure a climate signal; 
including economic development, the expansion of irrigation, 
wider implementation of modern cultivars, technological 
improvements, shifting international trade (Schmitz et al., 
2012), and the increasing and potentially changing demands of 
a growing population. Climate change will nevertheless be an 
additional stress on a country already facing considerable natural 
and societal challenges to meet the needs of a large and growing 
population. National production change projections obscure the 
more acute sub-regional impacts shown in this study that are 
likely to dominate local production trends unless targeted 
adaptations are implemented. 

Simulations extending out to the 2080s period (conducted as 
part of Yu et al., 2010) also shed light on the long-term 
trajectory of climate change impacts in Bangladesh. River floods 
were not modeled for this period, but temperature, precipita- 
tion, and sea level rise continue to have an increasingly negative 
effect on agriculture while the beneficial effects of C0 2 begin to 
level off. Additionally, the yield response to higher temperatures 
is not linear; temperatures exceeding critical thresholds during 
key phenological stages can trigger large production drops or 
even crop failure. The acceleration in detrimental temperature 
impacts on wheat is particularly great between the 2050s and 
2080s period. Thus, even areas in which climate change is 
projected to cause increased production in the early 21st 
century may face rapidly declining production due to climate 
factors in later decades if adaptation is not sufficient. Continuing 
research in agricultural observations, crop model development, 
and climate model projections will improve our ability to 
accurately assess these critical processes. 

Differences between models and emissions scenarios are 
substantial, but overall clustering of results suggests that the 
projections are robust and not strongly prone to divergent 
sensitivities. Additional uncertainty is introduced through ele- 
ments that could not be explicitly included in the integrated 
biophysical factor modeling framework, including: 

1. the effects of data-filling processes for climate scenarios; 

2. the assumption of stationarity in the statistics of extreme 
climate events; 

3. the use of generalized, high input agricultural management 
characteristics across Bangladesh; 

4. the sensitivity of production to C0 2 fertilization, high tempera- 
ture thresholds, and limited irrigation in the CERES crop models; 
and 

5. the neglected salinity effects and limitations in the simulated 
coastal domain. 


The results of this study therefore should be interpreted 
through the relative signs and magnitudes of changes for different 
crops, time periods, and climate scenarios. The spread between 
GCM-based climate scenarios and emissions scenarios has also 
been shown to inform the relative uncertainty of impacts when 
conducting risk assessments. 

Time, resources, and data availability limited the full imple- 
mentation of other considerations that would dramatically 
improve the quality of this (or future) multi-factor climate 
change impacts assessment. These include updated sea level rise 
scenarios including one with more rapid ice melt (e.g., Horton 
et al., 201 1 ), more realistic irrigation constrained by groundwater 
availability (Shahid, 2011), routines for salinity damage from 
irrigation and saltwater intrusion, the use of a weather generator 
to synthetically increase iterations and better capture extreme 
events, and climate change interactions with day-to-day and 
interannual climate variability. Continuing biophysical model 
development and data collection in Bangladesh will improve 
future analyses. 

The modeling framework used in this work allows for testing 
interactive adaptations, helping to quantify the “benefit” in cost- 
benefit prioritization. For example, the effects of virtual crops with 
alternative genetic attributes increasing climate resilience, the 
construction of additional polders, and coastal levies may all be 
simulated in this biophysical framework to gauge the impact on 
national production. 

The impact of climate change on the agricultural sector of 
Bangladesh is also highly dependent on climate impacts on food 
production in surrounding countries and the major agricultural 
areas that determine commodity prices. For example, projected 
impacts in this study have a different implication for Bangladesh if 
a large portion of other rice-producing countries experience 
production increases from regional climate changes, and likewise 
Bangladesh’s market share could increase if other rice-producing 
regions are severely impacted. Understanding such effects is one of 
the objectives of the Agricultural Model Intercomparison and 
Improvement Project (AgMIP; Rosenzweig et al., in press), which is 
now underway and includes a focus on South Asia. 
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